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1. Introduction 

Deep-inelastic electron-proton scattering experiments at HERA are measuring the structure 
function F2{x,Q'^) in the previously unexplored small x regime, 10~^ < a; < 10^'^. As usual 

= —q^ and x = Q'^/'^p ■ q, where p and q are the four momenta of the incoming proton and 
the virtual photon probe respectively. The deep-inelastic experimental observables reflect the 
small X behaviour of the gluon, which is by far the dominant parton in this kinematic region. 
In particular the small x behaviour of F2 is driven by the g ^ qq transition. At such small 
values of x soft gluon emission and the associated virtual gluon corrections give rise to powers 
of aslog{l/x) which clearly have to be resummed. Technically, in the physical gauge, this is 
equivalent to the summation of gluon ladder diagrams (with the virtual contributions leading to 
so-called gluon Reggeisation [1-6] or, alternatively, to the introduction of the non-Sudakov form 
factor [7-9]). These ladder diagrams are a universal feature of all small x processes driven by the 
gluon. For instance they occur in the perturbative QCD description of the structure functions 
F2 and Fl, heavy quark-pair and J/ip production, prompt photon production, deep-inelastic 
diffraction and deep-inelastic events containing a measured jet. 

The resummation of the leading powers of Q;slog(l/x) leads to the Lipatov (or BFKL) 
equation [1-10]. Here we wish to emphasize two characteristic features of the gluon obtained on 
solving this equation. First the gluon exhibits an x^^ growth as x decreases with A ~ 0.5. This 
power behaviour of the gluon should manifest itself in all the small x processes listed above 
since they incorporate the universal gluon ladder. We speak of a "hard QCD" or "Lipatov" 
pomeron with intercept 0^(0) = 1 + A ~ 1.5, applicable to values sufficiently large to be in 
the perturbative regime. This is to be contrasted with the universal "soft" pomeron (associated 
with ~ 0) whose intercept has been phenomenologically determined to be ap{0) ~ 1.08 |1TT|1 . 
The second feature concerns the transverse momenta, kx, of the gluons along the ladder. It is 
no longer true that the transverse momenta are strongly ordered as is the case in the "large" 
X, large "Altarelli-Parisi" regime. Rather, as we evolve down in x, we have diffusion arising 
from a "random walk" in which leads to a broadening of the starting fcy-distribution both 
to larger and to smaller values of k^- 

Instead of simply concentrating on the gluon distribution alone, here we study the impact 
of these perturbative QCD effects on the small x behaviour of the structure function F2{x, Q"^). 
However the implementation of this programme requires knowledge of the gluon for all k^ 
including the infrared region with k^ close to zero. In this intrinsic non-perturbative or con- 
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finement region the BFKL equation, which is based on perturbative QCD, is not expected to be 
vahd, and has to be modified. The major aim of this paper is to review and to critically examine 
the solutions of the Lipatov equation and their impact on the structure function predictions 
using physically motivated modifications for small k^. In particular we study the sensitivity 
of the predictions to a detailed parametrization of the infrared region which satisfies the gauge 
invariance constraints as — > 0. We find that the magnitude of F2 has a sizeable uncertainty 
whereas the x~'^ shape remains remarkably stable, with A being only weakly dependent on the 
choice of values of the infrared parameters. A comparison of the measurements of F2 with the 
predicted behaviour therefore provides a test of the Lipatov perturbative QCD effect. The 
Lipatov pomeron may, however, only give part of the increase of F2 and so we must allow for 
the effect of the "background" contribution when comparing with the data. 

The x~'^ growth of the gluon will eventually be suppressed by shadowing effects. We shall 
find these effects are relatively unimportant in the HERA regime, unless very small gluon 
concentrations exist within the proton. 

The outline of the paper is as follows. In section 2 we review the present knowledge of 
the properties of the solutions to the BFKL or Lipatov equation. In particular we illustrate 
the effects of diffusion and of the introduction of cut-offs in kT- In section 3 we introduce our 
infrared modifications of the BFKL equation and give numerical predictions for the structure 
functions. We confront these predictions with recent HERA data. Section 4 contains our 
conclusions. 



2. The behaviour of the solutions of the BFKL equation 

A central problem in small x physics is the stability of the solutions of the BFKL or Lipa- 
tov equation [1-9] to contributions from the infrared and ultraviolet regions of the transverse 
momenta of the emitted gluons. This is refiected in the dependence of the solutions to the 
choice of the transverse momentum cut-offs. Several general properties of the solutions of the 
Lipatov equation are known, which are scattered widely in the literature p|, |], ||, 0, |12|, [13|. 



In this section we draw these together and attempt to present a reasonably self-contained and 
coherent discussion. 

(a) Solution for fixed a.^ with no kx cut-offs 

It is helpful to begin by recalling the solution of the BFKL or Lipatov equation in the fixed 
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coupling case in the absence of cut-offs. The equation may be written in the form 



-X- 



d f{x,k') 
dx 



3a 



°° dk 
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f{x,k'')-f{x,e) , f{x,e) 



\k'^ - k^ 



+ 



{4:k'* + k*y- 



= K®f, (1) 

where /(x, A;^), the unintegrated gluon distribution, gives the probabihty of finding a gluon in 
the parent hadron with longitudinal momentum fraction x and transverse momentum squared 
k"^. To be precise, / is related to the more familiar integrated gluon distribution, g{x, Q^), by 

rQ' dk^ 



The BFKL equation sums up the leading powers of log(l/a;), that is it corresponds to the 
LL{l/x) approximation. In the genuine LL{l/x) approximation the strong coupling as should 
be taken as a fixed parameter. Although the integration over the transverse momentum squared 
in (1) does not contain any cut-off parameters, it should be emphasised that it is free from both 
infrared and ultraviolet divergences. However the solution contains infrared and ultraviolet 
singularities which will manifest themselves as non-trivial anomalous dimension(s); we amplify 
this comment in the discussion below (19). 

Since the kernel K is scale invariant, the equation can be diagonalized by the Mellin trans- 
form 



C+tCXD 



ZTTZ J c—ioo 



f{x,uj) 



{k'rf{x,u;)du; (2) 

{k^)-^-^f{x,k'')dk'' (3) 

where the integration contour in (2) is specified below. If we substitute (2) into (1) and change 
the variable k'^ to k'^u, then we obtain 

df{x, u) 



— X 



where 



K{u;) 



dx 

n Jo u 
3a 



= K{uj)f{x,u) 



(4) 



u"" - 1 
\u-l 



+ 



(41^2 + 1) 2 _ 



TT 



[2^(1) - ^{u) - ^(1 - a;)]. 



(5) 



where ^ is the logarithmic derivative of the Euler gamma function: "^{z) = r'{z)/r{z). If we 
substitute the solution of (4) into (2) we obtain 



fi^,k') = — / du;{k'rf{xo,uj) 



Xo 



(6) 



where xq is chosen small enough to ensure the validity of the Lipatov equation and yet large 
enough to see the effects of the evolution as x decreases from Xq. The integration contour in 
formula (6) should be located between the left-hand singularities of K{uj) arising from \l/(a;), 
and the right-hand singularities arising from ^'(l — ui) so that (6) reproduces the iterative 
solution of the original equation (1). It is therefore convenient to choose c — ^. Indeed K{u!) 
is symmetric about oj — ^ and, as it varies along the contour [uj — with — oo <v< oo), 

has its maximum value at this point. Thus, from (6), we see that the region OJ \ dominates 
the small x behaviour of f{x, k'^). We therefore write u — and expand K about v — 

k{\ + iv) = A-|AV + 0(i/^) (7) 

with 

A = ^41og2 (8) 

TT 

A" = — 28C(3) (9) 

TT 

where the Riemann zeta function ^(3) = 1.202. If we also expand the various terms in the 
integrand of (6) about u = ^ (that is = 0), then we find the behaviour of the gluon for 
X <^ Xq is given by the Gaussian integral 

fix,e) . ©"V>o4)^/_>-xp(-|A''log(^).^ + ..log(f;)) 
\xq) [27rA"log(a;o/a;)]^ ^ \2\"\og{xQ/ x) ) 

where 



logik") = i — {\og}{xo,\ + iv)%=o. (li; 



d 

du ' 

We have reproduced the solution of the gluon distribution originally obtained by Lipatov et al. 
with its characteristic x~^ behaviour with A given by (8), modulated by a (log(l/x))~^ factor. 

Formula (10) also displays exphcitly the diffusion pattern of the solution of the BFKL 
equation, that is a Gaussian distribution in log(/c^) with a width which grows as (log(l/x))5 
as X decreases. The position of the maximum of the Gaussian distribution (given by log(A;^) 
of (11)), as well as the normalisation of the solution, is controlled by the boundary conditions, 
that is by /(xq, k"^)- The rate of diffusion, however, is independent of the boundary conditions. 

The approximate analytic solution (10) only applies for x -C Xq. If we were to include the 
exp(— iAi/^) contribution in the expansion of f{xo, cu) about i/ — 0, then the A"log(xo/a;) factors 



in (10) would become A"log(xo/x) + A. We assume, simply for the purposes of illustration, 
that this modified form applies for all x < Xq. Fig. 1 shows the width of the Gaussian log 
/c^ distribution of f{x,k'^)/{k'^)^ as x decreases from xq. At Xq the "width" is given by the 
boundary conditions /(xq, k'^)/{k'^)^, though in practice this input distribution will not have a 
perfect Gaussian form in log fc^. In Fig. 1 we use dashed curves to emphasize the approximate 
nature of the treatment for x ^ Xq. It should be noted also that in a realistic treatment 
we find that the gluon distribution /(a;, A;^) samples k"^ uncomfortably close to the infrared 
(non-perturbative) region. 

(b) Physical examples of diffusion of the BFKL solution in fc^ 

The diffusion in log k"^ with decreasing x is a major problem in the applicability of the 
BFKL equation since it can lead to an increasingly large contribution from the infrared and 
ultraviolet regions of k"^ where the equation is not expected to be valid. We may illustrate 
diffusion using two physical examples from deep-inelastic electron-proton scattering. 

Given the unintegrated gluon distribution /(x, fc^) we can, in principle, calculate the be- 
haviour of the deep-inelastic structure functions -F2,l(x, Q"^) at small x through the so-called /ct 



factorization theorem |T^, |T5|. Then 

F<M') = /f^/'|^/(|;.*«)^■f(x',^='^Q=) (12) 

with i = 2, L. Symbolically we may write F = / ® F^°\ see Fig. 2(a), where / describes the 
gluon ladder and F^^^ the quark-box amplitude for gluon-virtual photon fusion. It should be 
noted that the integration over A;'^ extends down to k'"^ = and so knowledge of f{x/x', k'"^) in 
this region is, in principle, necessary for getting absolute predictions for -F2^l(x, Q"^). We return 
to this important point in Section 3. 

To illustrate the effect of diffusion in k"^ we use the LL{l/x) approximation to simplify (12) 

to 

F,{x,Q') = J —f{x,e)B,{k^Q') (13) 
where the "impact factors" Bi are 

B,ie,Q') = f%Fi'\x\k\Q% (14) 
Jo x' 

Now we may equally well rewrite the convolution (13) by factorizing at an intermediate link x\ 
along the gluon chain in Fig. 2(a) 

F.M^) = f^f{x„e)fJ^,kA (15) 



^4 ^ V--- Vxi 



where is a solution of the BFKL equation but with the boundary condition fixed at the 
"upper" end of the chain by the quark-box impact factor Bi{k'^,Q'^). The diffusion pattern is 



now determined by boundary conditions at both ends of the gluon ladder |T2|. To be specific, 
it is given by 

/ y^log(xo/xi)log(xi/x) 

( _ \og\e/ki) \ 

^ \ 2A"log(a;o/a;i) 2A"log(a;i/x) ) ^ ' 

where is determined by Bi{k'^,Q'^) and so k"^ ~ Q^. The variation of the width of the 
diffusion pattern, as xi varies between x and xq, is sketched in Fig. 2(a). Even for large Q^, the 
boundary conditions at Xq mean that the infrared region is penetrated leading to uncertainty 
in the predictions for Fi{x, Q^). 

This problem is overcome for deep-inelastic (x, Q"^) events containing an energetic measured 
jet (xj,/c|), see Fig. 2(b) [16-20]. We then have a 5{k'^ - fc|) distribution at the "bottom" of 
the gluon ladder and k'j can be chosen sufficiently large such that /(x, fc^) does not diffuse 
appreciably into the infrared region for physically accessible values of x/xj, see Fig. 2(b). 

(c) The cut-ofF dependence of the BFKL solutions 

One way to circumvent the diffusion into the non-perturbative region of k'^ is to (artificially) 
introduce an infrared cut-off fcg in the integration in the BFKL equation (1), and to explore the 
sensitivity of the solutions to reasonable variations of the value of the cut-off. In this subsection 
we show that if we introduce an infrared (or an ultraviolet) cut-off then the structure of the 
solution receives only minor changes, although the normalisation is altered. Most importantly 



the leading small x behaviour {x/xq)~ remains intact JTS . 

To investigate the sensitivity of the small x behaviour to the choice of cut-off, it is useful to 
transform /(x, k"^) to moment space 

h{n,k^) = [' dzz''-^f{x,k^) (17) 
Jo 

where z = x/xq. Then the small x behaviour is controlled by the leading singularity of h{n, k^) 
in the complex n plane. First let us recover the structure of solution (10) in the absence of a cut- 
off. In this case the equation, after its diagonalisation, (4) , by the Mellin transform, reduces to 
a simple algebraic equation in moment space. To be precise, we transform the integral version 
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of (4) to moment space and we find an algebraic equation witfi solution 



h(n,.) = 'f% (18) 
n — 1 — K[ijj) 



where fo{u}) = f{xQ,u!), and where h is related to h of (17) just as in (2) 

h{n,k^) ^ — r {k^Yh{n,u)duj. (19) 

It follows from (19) that the leading behaviour of h{n, k"^) as A;^ — > (/c^ — > oo) is controlled by 
the nearest singularity u!+{n) (a;_(n)) which lies to the right (left) of the contour of integration 
in the cu-plane, that is 

h{n, k^) - {k^f+ as k^ ^0 (20) 

(and as (A;^)'^- as A;^ — > oo). These u;± singularities come from the zeros of the denominator in 
(18), and the values u!± are equal, by definition, to the anomalous dimensions. Now, as we shall 
see in section 3, gauge invariance requires the driving term in the Lipatov equation to behave 
as 

fixo,k^) - k'^ as k'^ 0, 

corresponding to a pole at a; = 1 in /o(<^). In other words the presence of the anomalous 
dimension a;+ changes the small k"^ behaviour of h{n, k'^) from k"^ to (A;^)'^+(") where 0J+{n) < 1. 
It is in this way that the infrared singularities of the Lipatov equation manifest themselves. 
To determine the values of 0J±{n) we recall from (7) that in the neighbourhood of u — \ 

k{u) ^ A + iA>-i)^. 

Thus it follows from (18) that h{n,uj) has two nearby poles at 



uj = \±^J2{n-l- X)/\" = cu±(n). (21) 

These poles move together and pinch the contour in (19) when n = 1 + A and hence lead to a 
singularity in /i(n, /c^) at this point. Completing the contour in the left- or right-hand a;-plane, 
according as k"^ is large or small, gives 



^ ^ uj^-uj^ 2Vn - 1 - A ^ ' 
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If we insert this result in the inverse relation to (17), fold back the contour in the n plane 
to circle the square root branch cut and carry out the n integration, then we again find the 
behaviour of the gluon distribution f{x, k"^) given by (10); namely 



/ - ^-^[log(l/;^)]-^ (23) 

for small z = x/xq. 

To apply this technique to the solution of the Lipatov equation with a k"^ cut-off we express 
h{n,uj) as the sum of two components 

h+(n,uj) = deiey'^-^hin^e), (24) 

h-{n,uj) = / de{k^)-'^-^h{n,k'^), (25) 

cf. (3). The component functions h'^{h~) therefore have singularities in the right- (left-)hand uj 
plane respectively. In the case of an infrared cut-off, the algebraic equation which results from 
taking the moments of the integral version of (4) is 

h^{n,u) + h ln,uj) = -H ; (26) 

n — 1 n — 1 

and so 

h-M = ^'"'-'"-^f '"•"> . (27) 

n — 1 — K[uj) 

as compared to (18). In the cut-off case we have an additional constraint. We must adjust 
h^{n,uj) so that h~{n,u}) is free from singularities (i.e. poles) in the right half uu plane. In 
particular the numerator of (27) must contain a factor 



uj-uj+{n) = uj- \ - ^2{n-l-\)/\" (28) 
to cancel the pole at a; = ojj^in). As a consequence the pinch of the contour of integration 



in (19) no longer occurs. We see that h{n,uj) (and hence h{n,k^)) contains a y/n — 1 — \ 
singularity, rather than the l/\/n — 1 — \ singularity of the non cut-off case. Now when we 
fold the contour in the n plane round the branch cut and carry out the n integration we find 
the behaviour 

/ ~ z-^[log(l/^)]-i (29) 
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for small z = x/xq- (Clearly we obtain the same behaviour if we have an ultraviolet, and no 
infrared, cut-off - the roles of and are simply interchanged.) Finally we note that the 
infrared cut-off eliminates infrared singularities and so hin^ k^) ~ A;^ as A;^ — > rather than the 
anomalous behaviour shown in (20). 

So far the leading small x behaviour, x~'^, has remained intact to the multilations of the 
Lipatov equation. However if we introduce both an infrared and ultraviolet cut-off then the 
exponent A becomes cut-off dependent. 

(d) Fixed and running a^: the eigenvalue spectrum of the BFKL kernel 

The exponent A controlling the small x behaviour, a;"'^, of the gluon is given by the maximal 
eigenvalue of the BFKL kernel, K{uj); see (6-8). If we use a fixed value of as in the BFKL 
equation, (1), then the eigenvalue spectrum is continuous, and remains so in the presence of 
either an infrared or ultraviolet cut-off. Moreover the maximum eigenvalue (the branch point 
of the cut) does not change. 

The situation is different if we introduce both an infrared cut-off fcg and an ultraviolet cut off 



^max F3| ■ Then the eigenvalue spectrum becomes discrete. The maximum eigenvalue, Amax, and 
the separation between the eigenvalues can be shown to depend on the quantity t = /Cmax/^o- 
For large log t the distance between the eigenvalues becomes proportional to 1/log t and 

Amax = A,nax(t = 00)+0(l/l0gt), (30) 

SO in the limit t — oo we do indeed recover the continuous spectrum. 

So far we have considered a fixed value of the coupling, as, in the BFKL equation, (1). 
There are indications that we should allow the coupling to run, that is as — > asik"^)- Here we 
make this physically reasonable replacement, although, as yet, there is no rigorous proof. At 
present the main argument is that in the strongly ordered fc^ limit the BFKL equation will 
produce the Altarelli-Parisi equation in the double leading logarithm approximation if we take 
as = as(fc^). The introduction of running has the effect of suppressing the importance of 
the ultraviolet cut-off and enhancing the dependence on the infrared behaviour. Moreover in 
this case, even with no ultraviolet cut-off, the eigenvalue spectrum of K^u) is discrete, with 
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Amax sensitive to the choice of the infrared cut-off k^. 



3. Numerical solution of the BFKL equation 

(a) Treatment of the infrared region 

We shall use a running coupling asik"^) in the BFKL equation and so we will need to focus 



attention on how to deal with the infrared region. The simplest procedure is to introduce 
a cut-off /cq (as in section 2(b)) so the BFKL equation becomes 



df 3a, (P) f^^dk'^ 



X-— = k 



dx 71 Jkl k'"^ 



f{x,k'')-f{x,e) , f{x,e) 



+ 



(31) 



where (for simplicity) the same cut-off is used in the real emission term and in the virtual 
corrections. To calculate F2 we would impose the same cut-off on the convolution integral (12) 
which occurs in the fc-r-factorization theorem. 

The above cut-off which completely eliminates the infrared region k^ < k^ is rather drastic. 
Clearly a better procedure which incorporates this region, at least in an approximate way, is 
desirable. The problem is that the BFKL equation is not expected to be valid when the gluon 
momenta enter the non-perturbative region of small k"^. One way to overcome the problem 
is to introduce non-perturbative (albeit phenomenological) gluon propagators which are finite 
at = p2| , p3| and hence to eliminate the potential infrared singularities of the solution. 
Alternatively we could modify the BFKL equation by explicitly subtracting the point fc^ = 
]15| . Or, as we shall do below, we could explicitly introduce a factor 1 — G{k'^) into the 



solution with G(0) = 1. The common feature of these modifications of the infrared region is 
that the solution of the Lipatov equation (and the driving term) vanish as A;^ as /c^ — 0. This 
requirement that 

/(x,A;^) ~ as fc^ ^ (32) 

is a consequence of gauge invariance or to be precise of the colour neutrality of the probed 
proton [^, p!0[| . 

In the present paper, we attempt to go beyond the simple k'^ cutoff approximation and to 
model the low k'^ region in a more systematic fashion. We assume that the small-A;^ behaviour 
of the gluon distribution is driven by a form factor G{k'^) such that 

f{x,k^) ~ Const. [1 - G(fc2)] (33) 
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for k"^ 0. We take 

where the parameter is related to the radius of the gluonic form factor of the proton. If this 
is taken to be of the same magnitude as the radius characterising the hadronic electromagnetic 
form factor then we would have A;^ ~ 0.5 GeV^; however estimates based on the QCD sum rules 
prefer a larger value, fc^ ~ 1 — 2 GeV^. 

We then proceed by splitting the integration region for real gluon emission (the term in- 
volving f{x, k'^)) in (31) up into two parts, namely 

region (A): to /cq 
region (B): /cq to oo. 

In region (B) the BFKL equation as it stands is taken to hold. In region (A) we assume that /cq 
is sufficiently small that the behaviour given in (33) is a good approximation. If we parametrise 
f{x, k'"^ < fcg) in this form, then the integral (A) can be calculated analytically, and in this 
way we have a physically motivated approximation for the infrared contribution to the BFKL 
equation (31). As a further modification to the low-fc^ region we 'freeze' the argument of as by 
using as{k'^ + a^), with = 1 GeV^, in both the evolution equation and in the factorisation 
formula (12) used to calculate F2 and F^. 

The modified BFKL equation can then be used to evolve the gluon distribution down in 
X starting from a suitable input distribution f{xQ,k'^). This boundary condition must be 
consistent with Altarelli-Parisi evolution for large fc^, that is 



PAP^ ^2^ d{xog{xo,Q'^)) 



9iogg2 



(35) 

Q2=fe2 



where we take g{xo, Q ) from the Altarelli-Parisi evolution of MRS partons [24|. The boundary 



condition must also have a small k'^ behaviour consistent with our approximations, so we take 

f{xo,e) = {l-G{e))f^^{xo,e + a') (36) 

for k^ > /eq, whereas for k^ < k^ we "freeze" the evolution of /^^ at k^ + a^. The parameter 
a} = \ GeV^ just softens the low-/c^ behaviour of f^^ which tends to be unreliable. The 
important fact is that the input distributions approach /^^(xq, A;^) for large fc^, as they must, 
and also embody a suitable behaviour for smaller k^ . 
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In preparation to see the diffusion in A;^ develop as we proceed to small x, we plot the 
boundary conditions in the form f{xQ,k'^)/{k'^)^ as suggested by (10). Sample distributions 
are shown in Fig. 3 for different choices of k"^ and k^. We see the input distributions have an 
approximate Gaussian form in log k"^. Fig. 4 shows the evolution of the distribution, for the 
choice kl = kQ = 1 GeV^, as we proceed to smaller x using (31) and (33). We see both the 
diffusion to large k"^ and the x~^ type growth. There is no diffusion into the infrared region 
since we impose the phenomenological form (33). The diffusion to large k'^ is more apparent in 
Fig. 5 which shows the distributions of Fig. 4 normalised to a common value. Even in the limit 
of very small x, the rate of this diffusion will differ from that given by the analytic form (10), 
since here we are using a running coupling a^- 

(b) Consistency constraint on k^ 

The numerical solution f{x, k"^) of the Lipatov equation, modified as above, is found to be 
much more sensitive to the choice of k^ than to /cg, see ref. [25]. Therefore here we choose 
/cq = 1 GeV^, and we concentrate on investigating the sensitivity of the results to variations 
of k"^. However as may be anticipated from the discussion in section 2 it is the magnitude of 
/, and not the shape, which is particularly sensitive to fc^. That is, with decreasing x, an 
behaviour sets in with a numerical value of A only weakly dependent on fc^ |25|. 

There is a consistency requirement between the "input" and "output gluon" which, in 
principle, can be used to estimate the value of kl- The constraint is that the gluon distribution 
calculated from 

xg{x,Q') = ^f{x,k'), (37) 

the inverse relation to (35), should match the phenomenological input gluon distribution. The 
comparison is shown in Fig. 6 for different choices of k'^. Here we use as input at xq = 0.01 a 



gluon satisfying the leading order Altarelli-Parisi equations (the dash-dotted curve), but 



based on the Dq type parametrizations of ref. |26|. Fig. 6 shows that there is good agreement 



between the "input" and "output" gluons for k"^ ^ 1 GeV^. (In Fig. 6 we also compare our 
input gluon with the gluons of the Dq and D'_ next-to-leading order analyses of ref. ||2^). In 
practice, it would be misleading to impose this constraint too rigorously. The input gluon is 
not well known at Xq = 0.01, particularly at the low values of which are necessary for this 
comparison. Nevertheless it is encouraging that the estimate of k"^ appears reasonable, and 
suggests that we should consider QCD predictions with k^ chosen in a range of about 0.5 to 2 
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(c) Numerical predictions for F2 and at small x 

We may use the solution f{x, k"^) of the modified Lipatov or BFKL equation to predict the 
structure functions Q^) at small x via the factorization theorem (12). The factorization 
formula has the symbolic form 

i^Lip = f^F^ (38) 

and gives the contribution to Fi arising from the BFKL resummation of soft gluons. Recall 
that describes the quark box (and "crossed" box) photon-gluon fusion process shown in 
the upper part of Fig. 2(a). The gluon is off-mass-shell with virtuality approximately equal to 
k'"^. The explicit formula for Ff\x\ k'"^, Q^) can be found in refs. ||l^ 

Before we can obtain a realistic estimate for the structure functions Fi we must include 
the background non-Lipatov contributions F^^. A reasonable choice at small x is to assume 
that F^^ gradually increases like x~^'^^ with decreasing x (as might be expected from a "soft" 
pomeron with intercept ap(0) = 1.08 ]ll|). To be precise we take 

F^^{x,Q') = Fj'^{xo,Q'){x/xo)-'-'' (39) 

with xq = 0.1. The values we use 

fori^^s(xo,Q2) 

are listed in the figure captions. The resulting 
predictions for the small x behaviour of Fi = F^"^^ + fJ^^ with i = 2,L are shown in Figs. 7-10 
for various values of Q^. In each figure the continuous curves show the predictions for two 
choices of the infrared parameter k"^, namely k"^ = 1 and 2 GeV^. 

The recent HERA measurements of F2 are shown on Figs. 7 and 8. There is general 
agreement between QCD and the data. In particular they both show a dramatic increase with 
decreasing x, and lie well above a straightforward extrapolation of the fixed target measurements 
that exist for x > 10~^ pi], |33]. Certainly the data indicate support for the x~^ type behaviour 



arising from the BFKL leading log(l/x) resummation. But before we draw conclusions we must 
consider the effects of shadowing corrections. 

(d) Inclusion of shadowing 

The growth of the gluon density with decreasing x means that there is an increased proba- 
bility that the gluons will interact and recombine. To allow for the effects of this recombination 
or parton shadowing we incorporate an additional term in (31) 

_,dfix^ = K^f-j^yAk')[xg{x,k-)r (40) 
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where g is given by (37). The additional term, quadratic in g, in (40) is the leading order 
shadowing approximation; the negative sign leading to a suppression in the growth of the gluon 
density with decreasing x. The crucial parameter is R, where nR^ specifies the transverse area 
in which the gluons are concentrated within the proton. 

For illustration we take R = 5 GeV~^ (corresponding to gluons uniformly spread across 
the proton) and R = 2 GeV~^ (assuming the gluons are concentrated in "hot-spots" within 
the proton). The dashed curves in Figs. 7-10 show the effect of these two shadowing scenarios 
respectively on the = 1 GeV^ prediction. We note that the shadowing effects are now 
slightly stronger than in the case when the region of small k'"^ was entirely neglected. 



However shadowing is still rather a weak effect in the HERA regime and sets in very gradually, 
unless compact "hot-spots" of gluons occur. In particular we are far from the saturation limit. 
We have found that F^^^ behaves like Cx~^ for x < 10^^ where the predicted value of A 



is relatively insensitive to the uncertainties associated with the infrared region. The inclusion 
of shadowing means that A is no longer constant but that its value decreases with decreasing x, 
as illustrated by the dashed curves. The predictions for = 15 GeV^ (and k^ = 1 GeV^) are 
shown in Fig. 11. Conventional shadowing {R = 5 GeV^^) has relatively little impact on A, and 
even if "hot-spots" were to exist A remains significantly above the soft-pomeron expectation of 
0.08. 

We have checked that the large values of A, which are essentially independent of x for x < 
10~^, remain true for other physically reasonable choices of the background F^^{x, Q"^). Indeed 
the values of A obtained in Fig. 11 are much larger than would result from any straightforward 
extrapolation of the fixed target F2 data to small x. Consider for example the value of A 
obtained from the extrapolation of the D'^ set of MRS partons ||2^ to small x. (The Dq partons 



provide an excellent description of the fixed target data which only exist for x > 10 ^). The 
parton set Dq develops, via Altarelli-Parisi evolution, a small x behaviour of the form 

F2~exp[2(e(g^,g^)log(l/a;))^/2] (41) 

where the "evolution length" 

5,gg,«'^).rf!«!!-W!) ,42) 



Ql Q 



TT 



The increase of F2 with decreasing x can be translated into an effective x ^ behaviour in the 
HERA regime. In fact the value of A is found to be about 0.11 at = 15 GeV^ and 0.15 at 
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= 30 GeV^; slightly increasing with on account of the increase in the evolution length 
C(Qo' Q"^)- Note that these small values of A rely on the choice of a sufficiently large value of 
Ql, for instance MRS evolve from Ql = A GeV^. 

However GRV |Q have obtained partons by evolving from a valence like input at Qg = 0.3 
GeV^. The very low value of Ql corresponds to a relatively large evolution length ^{Ql, Q^) for 

in the HERA region i.e. ~ 20 GeV^. In this way they obtain a steeper small x behaviour 
for F2 (see (41) and (42)) compatible with the data. In fact this double leading log. behaviour 
mimics an form with A ~ 0.4 in HERA regime. However we do not believe that this is an 
acceptable explanation of the data since the steepness is mainly generated in the very low 
region where perturbative QCD is invalid (see, for example, [0^). 

In order to obtain the steep x^^ type behaviour with A ~ 0.5 within the Altarelli-Parisi 
formalism (and when the evolution starts at moderately large value of Qo ~ 5 GeV^) one 
has to impose this steep behaviour in the parametrisation of the starting gluon and sea-quark 
distributions at the reference scale Qg as it is done for instance in the case of the D_ set of 
MRS partons [^, . In this procedure however one does not use the BFKL equation, that is 
the singular behaviour is not generated explicitly by QCD dynamics. 



4. Conclusions 

The recent measurements ^ of the deep inelastic structure function F2 at HERA 
explore the small x regime for the first time. The data show that F2 increases as x decreases 
from 10^^ to a few xlO~^, and do not follow a straightforward extrapolation of the fixed target 
measurements [^T], ^ that exist above x ~ 10^^. This novel behaviour is in line with the 
growth anticipated from perturbative QCD via the fcy-factorization formula, symbolically of 
the form 

F2 = f^Fi'\ (43) 

which links the small x behaviour of F2 with that of the universal unintegrated gluon dis- 
tribution / via the quark box contribution ^2°'' to photon-gluon fusion. The growth of F2 
with decreasing x is thus associated with the BFKL leading log(l/x) summation of soft gluon 
emissions which yields the small x behaviour /(x, k^)) ~ x^^ with A ~ 0.5. 

However the contribution from the infrared (low /c^) region, which occurs in the convolution 
of (43), leads to a sizeable uncertainty in the predictions. Motivated by the apparent agreement 
between the data and the expectations of perturbative QCD, we have attempted to improve the 
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treatment of the low regime. In particular, rather than imposing a low kr cut-off, we employ 
a physically motivated low kx form for f{x, k^) which allows us to extrapolate right down to 
k^ = 0. This reduces the uncertainty in the predictions for F2, although the normalisation 
still depends significantly on the choice of the value of an infrared (form factor) parameter 
kl- However the effective slope A which specifies the shape is much less sensitive to the 
ambiguities at low kx- In Section 2 we gathered together general arguments which suggested 
that the slope might be relatively immune from infrared effects and in Section 3 we performed 
explicit numerical tests to verify the result. 

In Figs. 7 and 8 we compared the perturbative QCD calculations for F2 with the recent 
HERA data. The dramatic growth with decreasing x is apparent in both the data and the 
QCD predictions. The various curves show the sensitivity of the QCD determination to the 
variation of the infrared cut-off and to shadowing effects. Fig. 11 translates these results into 
a comparison for A [where A is defined by the the growth of the BFKL component of 
F2{x, Q'^)]. We see that QCD indeed predicts an approximate behaviour at small x (or in 
the extreme case of "hot-spot" shadowing an type growth with decreasing x). 

The results from HERA are very encouraging and suggest that HI and ZEUS may have 
seen the first evidence for the BFKL growth arising from the leading log(l/a;) soft gluon re- 
summation. Figs. 7,8 and 11 should be regarded as an indication of what may be learnt when 
much higher statistics data become available and allow a detailed comparison. At the moment 
it is only possible to make predictions in the leading log(l/x) approximation. Much effort is 
being devoted to obtaining the next-to-leading contributions and, when available, these should 
be incorporated. Inspection of Fig. 11 suggests that, from a study of the x dependence of A 
for X < 10~^, we may then be able to quantify the effects of shadowing. 
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Figure Captions 

Fig. 1 The variation in the width of the Gaussian log fc^ distribution of /(x, as we 

evolve down in x below the starting value xq. 

Fig. 2 The upper diagrams show a gluon "ladder" contribution to small x (a) for deep-inelastic 
scattering and (b) for deep-inelastic scattering together with an energetic jet. The quark 
box factor F^'^^ implicitly includes the contribution of the crossed box. The lower sketches 
show the variation of the width of the log/c^ distributions of (16) as a function of Xi. 

Fig. 3 The boundary conditions for f /{k"^)^ at x = 0.01 used to solve the modified BFKL 
equation, for various choices of k"^ and fcg. {k^ = 1 GeV^ where it is not specified). 

Fig. 4 The evolution of f /{kP')^ as we step down in x using the modified BFKL equation with 
kl = kl = l GeV2. 

Fig. 5 As for Fig. 4 but with the distributions normalised to a common value at /c^ = 1 GeV^. 

Fig. 6 The self consistency of the gluon at x = 0.1. The dash-dotted curve is the input gluon 
2^ and the continuous curves show the output gluon obtained from f{x,k'^) via (37), 



where / itself is determined from the input gluon with different choices of k"^ (but with 
fcp = 1 GeV^). The dashed curves, which correspond to the Dq and D'_ gluons of ref. p7[] , 
are to illustrate the ambiguity in the input gluon. 

Fig. 7 The perturbative QCD predictions for F2{x,Q'^) at = 15 GeV^ obtained from the /ct 
factorisation formula, (12). The continuous curves correspond to the infrared parameter 
kl = 1 and 2 GeV^ respectively, with shadowing neglected. The dashed curves show 
the suppression caused by conventional (i? = 5 GeV^^) and "hot-spot" {R = 2 GeV~^) 
shadowing for the choice k^ = 1 GeV^. The data are from the HI p9| and ZEUS [0 
collaborations. The background contribution is given by (39) with Fp{xo) = 0.384. 

Fig. 8 As for Fig. 7 but for = 30 GeV^ with F^^ixo = 0.1) = 0.391. 

Fig. 9 The curves are as for Fig. 7 but for the longitudinal structure function Fl{x,Q'^ = 
15GeV^). The background contribution is given by (39) with F^^ = 0.04 at xq = 0.1 and 
g2 = 15 GeV2. 

Fig. 10 As for Fig. 9 but at = 39 GeV^ 
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11 The effective slope A, defined by F2 = ^2^^ + Cx~^ wfiere is given by (39), for various 
clioices of tlie infrared parameter k"^. Tlie lower two (dot-dashed) curves show the effect 
of the conventional (i? = 5 GeV~^) and "hot-spot" (R = 2 GcV^^) shadowing on the 
— 1 GeV^ predictions. The "data" points are calculated from the HI and ZEUS data 
shown in Figs. 7 and 8. 
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